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ABSTRACT 

We study the effect of cooling flows in galaxy clusters on the Sunyaev-Zel'dovich (SZ) distortion and 
the possible cosmological implications. The SZ effect, along with X-ray observations of clusters, is used 
to determine the Hubble constant, H . Blank sky surveys of SZ effect are being planned to constrain the 
geometry of the universe through cluster counts. It is also known that a significant fraction of clusters 
has cooling flows in them, which changes the pressure profile of intracluster gas. Since the SZ decrement 
depends essentially on the pressure profile, it is important to study possible changes in the determination 
of cosmological parameters in the presence of a cooling flow. We build several representative models of 
cooling flows and compare the results with the corresponding case of gas in hydrostatic equilibrium. We 
find that cooling flows can lead to an overestimation of the Hubble constant. Specifically, we find that 
for realistic models of cooling flow with mass deposition (varying m with radius), there is of the order 
~ 10% deviation in the estimated value of the Hubble constant (from that for gas without a cooling flow) 
even after excluding ~ 80% of the cooling flow region from the analysis. We also discuss the implications 
of using cluster counts from SZ observations to constrain other cosmological parameters, in the presence 
of clusters with cooling flows. 

Subject headings: Cosmology: cosmic microwave background, distance scale, large scale structure of the 
universe ; Galaxies: clusters: cooling flows 



1. INTRODUCTION 

Galaxy clusters are extensively observed in optical, X- 
ray and radio bands. In the radio band, a cluster can be 
observed in the the Rayleigh-Jeans side of the cosmic mi- 
crowave background spectrum, as a dip in the brightness 
temperature, due to Sunyaev-Zel'dovich effect (SZ) (Sun- 
yaev & Zel'dovich 1972; for a comprehensive review, see 
Birkinshaw 1999). The SZ distortion appears as a decre- 
ment for wavelengths > 1.44 mm (frequencies < 218 GHz) 
and as an increment for wavelengths < 1.44 mm. SZ effect 
has the advantage that the SZ intensity, unlike that of the 
X-ray, does not suffer from the (l + z)~ 4 cosmological dim- 
ming. As discussed by numerous authors (Birkinshaw & 
Hughes 1994; Silverberg et.al. 1997), one can combine the 
X-ray and radio observations for clusters to determine cos- 
mological parameters. This has been done in the recent 
years to determine the Hubble constant H (Birkinshaw 
1999). The SZ signal is, however, weak and difficult to de- 
tect. Recent high signal to noise detection have been made 
over a wide range in wavelengths using single dish obser- 
vations: at radio wavelengths (Herbig et al, 1995, Hughes 
& Birkinshaw 1998), millimeter wavelengths (Holzapfcl et 
al. 1997, Pointecouteau et al. 1999) and submillimeter 
wavelengths (Komatsu et al. 1999). Interferometric ob- 
servations have also been carried out to image the SZ ef- 
fect (Jones et al. 1993, Saunders et al. 1999, Reese et 
al. 1999, Grego et al. 2000). Other than estimating H 
combining SZ and X-ray data, SZ effect alone can also 



be used to determine the cosmological mass density f2 of 
the universe (Bartlett & Silk 1994, Oukbir & Blanchard 
1992, 1997, Blanchard & Bartlett 1998). However, these 
procedures generally assume the cluster gas is spherical, 
unclumped and isothermal. Almost all clusters, however, 
show departures from these simplistic assumptions with 
some to a large extent. 

Departures from these simple assumptions can lead to 
systematic errors in the determination of the different cos- 
mological parameters (Inagaki et al 1995), especially, it 
was seen that non-isothermality of the cluster can lead to 
a substantial error in values of the cosmological parame- 
ters. Temperature structure in a cluster can be the result 
of the shape of the gravitational potential (Navarro et al 
1997, Makino et al 1998), or it can arise due to the fact that 
the initial falling gas in the cluster potential is less shock 
heated than the later falling gas (Evrard 1990). In fact hy- 
drodynamical simulations of isolated clusters also show a 
definite temperature structure and can introduce error in 
the value of H when compared to the traditional isother- 
mal /3-models (Yoshikawa 1998). Roettiger et al (1997) 
have shown that cluster mergers can result in deviations 
from both sphericity and isothermality. Observationally, 
the main handicap arises from the fact that the thermal 
structure of clusters are hard to measure, and the temper- 
atures generally taken in analyses are the X-ray emission 
weighted temperature, usually measured over a few core 
radii. So, in general an isothermal description of the clus- 
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ter is taken (or sometimes a phenomenological tempera- 
ture model based on the Coma cluster: for example see 
Eq 73 of Birkinshaw 1999). 

In this paper, we study another important phenomenon 
that can substantially change the temperature structure, 
viz. a cooling flow. Cooling flows in clusters of galaxies 
(for an introduction, see Fabian et al 1984) is a well estab- 
lished fact by now, and it is seen that around 60 — 90% 
of clusters exhibit cooling flows in their core with 40% 
of them having cooling flows of more than IOOMq yr^ 1 
(Markevitch et al 1998, Peres et al 1998, Allen et al 1999). 
In the largest systems, the mass deposition rate can be as 
high as 1000M Q yr" 1 (Allen 2000). The idealised picture 
of a cooling flow is as following: Initially when the clus- 
ter forms, the infalling gas is heated from gravitational 
collapse. With time this gas cools slowly and a quasi hy- 
drostatic state emerges. However in the central region, 
where energy is lost due to radiation faster than elsewhere, 
an inward 'cooling flow' initially arises due to the pressure 
gradient (Fabian 1994). This can modify the SZ decrement 
and act as a systematic source of error in the determina- 
tion of the cosmological parameters. 

Schlickeiser (1991) has shown that free-free emission 
from cold gas in the cooling flow can actually lead to an 
apparent decrease of the SZ effect at the centre. Since 
the central cooling flow region is generally very small, the 
isothermal /3-model of cluster gas can still be used for the 
majority of cluster region even for cooling flow cluster, 
with the extra precaution of excluding the central X-ray 
spike from the X-ray fit, and a corresponding change made 
in the fitting of the SZ decrement. This is only possible for 
nearby clusters, however with well resolved cluster cores. 

Naively, the change in the central SZ decrement y(0) 
can be seen as follows : For a non cooling flow cluster, the 
central decrement is given by the line of sight integral of 
the electron pressure through the cluster centre along the 
full extent of the cluster. If the cluster has a maximum ra- 
dius r c i, then the central SZ decrement at RJ wavelengths 
can be written as y(0) = —^-^s Jn"' p e dl. For a cluster 
with a cooling flow, let us suppose that the electron pres- 
sure p ei drastically falls below a certain radius r s , which 
is typically well inside the core of the cluster. The re- 
sulting central decrement is then j/(0) ~ — 4^^ J r °' p e dl. 
Depending on the distance of r s from the cluster centre 
[r s Ri (0.1 to 0.3)r core ], there wiil be a change in the value 
of y(0) by w 5% — 25%. However, this simplistic view may 
not be true. This estimate assumes that the pressure pro- 
file remains a (3 profile outside the radius r s . The pressure 
profile, however, need not follow the (3 profile once cooling 
flow starts and it can deviate from it substantially even 
for radii much larger than r s . As a matter of fact, there 
can actually be an increase in the pressure for a large re- 
gion inside the cooling flow, before a sudden drop inside 
r s . Since the usual proecedure for estimating the Hubble 
constant depends on fitting (3 profiles to the SZ and X-ray 
profiles, to estimate r core , this change in the pressure pro- 
file due to cooling flow can distort the estimation of r core 
and hence, the value of H Q in a non-trivial way. We study 
this effect in detail in later sections. 

In this paper, we have investigated the problem of cool- 
ing flow induced change in the temperature and density 
profile, its effect on the SZ effect, and its subsequent effect 



on the determination of cosmological parameters. In §2 we 
briefly review the physics of Sunyaev Zel'dovich effect; §3 
is devoted to the physics of cooling flows and discussing the 
cooling flow solutions; in §4 we look at the effect of cooling 
flow solutions on SZ effect, both on the determination of 
H and Q ; we conclude in §5 with a brief comment on 
how this work differs from other work and the relevance of 
this paper. For the SZ effect, our notation and approach 
mainly follows that described in Barbosa et al(1996). 

2. DETERMINING HUBBLE CONSTANT WITH 
SUNYAEV-ZEL'DOVICH EFFECT 

2.1. The Sunyaev-Zel'dovich Effect 

The integral of the electron pressure along any line-of- 
sight through the cluster determines the magnitude of the 
distortion of the apparent brightness temperature of the 
cosmic microwave background (CMB) due to SZ effect. 
This is quantified in terms of the Compton y-parameter: 
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where k B is the Boltzman constant, T e is the gas tem- 
perature, m e is the electron rest mass, n e is the electron 
number density, c is the velocity of light and ot is the 
Thomson scattering cross-section. This occurs through 
the inverse-Compton scattering, by the hot intracluster 
gas, of the CMB photons propagating through the clus- 
ter medium, and the energy transfer in this interaction 
between hot electrons and CMB photons resulting in a 
distortion to the CMB spectrum. The SZ surface bright- 
ness at a position 9 of the cluster with respect to the mean 
CMB intensity is given by 



a: is a dimensionless frequency parameter 
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where h is the Planck constant, v is the observing frquency 
and To is the CMB temperature at the present epoch: 
T ~ 2.73K. The function j v {x) describes the spectral 
shape of the effect 
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Since the total photon number is conserved in the inverse 
Compton scattering process, upscattering of the photons, 
the spectral dependance gets an unique shape, through a 
decrement in the brightness temperature at lower frequen- 
cies while an increase is observed at higher frequencies. 

The Sunyaev-Zel'dovich effect provides a unique observ- 
ing approach to traditional methods — which use X-ray 
temperature and X-ray luminosity. The X-ray studies have 
a major disadvantage because of 'cosmological dimming', 
the surface brightness of distant X-ray sources falls off as 
(l + z)~ 4 , and for this reason, obtaining samples of clusters 
at cosmological distances is challenging. The SZ effect has 
the distinct advantage of being independent of the distance 
to the cluster. The SZ flux density from a cluster will di- 
minish with distance to the cluster as the square of the 
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angular-size distance; in contrast to X-ray flux densities 
from clusters which diminish as the square of the luminos- 
ity distance to the cluster. If the observations resolve clus- 
ters, particularly at lower rcdshifts, the observed sky SZ 
temperature distribution will be sensitive to the thermal 
electron temperature structure within the clusters; once 
again this may be contrasted with X-ray emission images 
of cluster gas distributions which are mainly sensitive to 
the gas density distribution. 



extension of the cluster. Conventionally, R c iuster — oo, 
and then from equations (5), (6), (10), (11), we get 
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2.2. Determination of Hubble constant 

The method for the determination of Hubble constant 
using Sunyaev-Zel'dovich effect uses two observable quan- 
tities : 1) AT/T of the CMB due to SZ effect; 2) the X-ray 
surface brightness Sx of the cluster. These can be written 

as 
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where r is the distance to the line of sight from the clus- 
ter centre, l ma x and l m in give the extension of the cluster 
along the line of sight, is the X-ray emissivity and dl 
the line element along the line of sight. The X-ray emis- 
sivity in the frequency band v = v\ to v 2 can be written 
as 

dLx 2 /m \ /rr\ 

n e a{T e ;v 1 ,u 2 ,z), (7) 
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In the above equations we have assumed primordial abun- 
dance of hydrogen and helium and have set X = 0.76, e 
is the electron charge, K — h/(2ir), u = 2-kUv /k B t e , and 
gff(T e ,u, Z) is the velocity averaged Gaunt factor for the 
ion of charge Zc (Kcllog 1975). Traditionally, to model the 
cluster gas distribution one takes the following density and 
temperature profiles (Cavaliere & Fusco-Femiano 1978) 
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T e (r) = T iso — constant, 



(10) 



(11) 



where n e o is the central electron density and r core the core 
radius of the cluster. The above expressions are used as 
an empirical fitting model, and the parameter '/?' is re- 
garded as the fitting parameter. The equation holds for 
< r < Rciuster, where Rciuster is the maximum 'effective' 
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where T(x) is the gamma function. Since both the central 
CMB decrement and the X-ray surface brightness are ob- 
served, one can then combine equations (12) and (13) to 
estimate the core radius as 
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In the above equation 8x,core is the angular core radius ob- 
served in X-ray, and T e ju is the X-ray flux averaged tem- 
perature (obtained from fitting the observed X-ray spec- 
trum to the theoretical spectrum expected from isothermal 
case). This X-ray emission weighted temperature is given 

by 

C"' r T e (r)an e 2 r 2 dr , s 
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The point to be noted is that r vir is the virial radius of 
the cluster and its choice depends on the observer. If the 
temperature has a spatial structure then the Ti SO inferred 
from such a procedure may give different values depending 
on how much of the cluster is taken in making the above 
average. It has been seen (Yoshikawa et al 1998), that this 
can lead to a substantial change in SZ effect inferred, and 
thus to the value H . 

The angular diameter distance dA can be approximated 
for nearby (z <C 1) clusters as 
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Thus finally we have got an estimate of dA{z) as 
r c ,est/6x,core- Now, if from other observations we know 
the cosmological parameters f2 and A, then one can esti- 
mate the Hubble constant 

As can be seen from the above equations, the value of H 
depends crucially on the many assumptions of isothermal- 
ity and /3-model density distribution of the cluster. Cool- 
ing flow changes both of these and so it can significantly 
affect the value of the Hubble constant. 
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3. COOLING FLOWS IN CLUSTERS 

3.1. Preliminaries 

From X-ray spectra of clusters, it is known that the 
continuum emission is thermal Bremsstrahlung in nature 
and originates from diffuse intracluster gas with densities 
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and temperatures around 10 — 10 K. 



The gas is usually believed to be in hydrostatic equilib- 
rium. If, however, the density in the inner region is large 
enough, so that the cooling time is less than the age of the 
cluster, then there is a 'cooling flow' (Fabian et al 1984 and 
references therein) . Of course there would be a flow only 
when the dynamical time is also shorter than the cooling 
time scale (t age > t coo i > t dyn ). 

The basic equations of cooling flow are: 
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where 9 — 2k bT '/ ' pm p , p, is the mean molecular weight, 
and m,p is the proton mass. For steady flows with constant 
mass flux, p = 0. This implies u = rh/4irpr 2 for steady 
flows. (Note that in cooling flows both u and m are neg- 
ative. The subscripts b refers to baryons and t refers to 
total i.e baryons + dark matter. However, we assume the 
baryonic contribution to the total mass negligible w.r.t to 
the dark matter contribution.) 

M{r) describes the distribution of the total mass and de- 
pends on the details of dark matter density profiles (see be- 
low). A(6) is the cooling function defined so that n e n p A(9) 
is the rate of cooling per unit volume. We use an analyt- 
ical fit to the optically thin cooling function as given by 
Sarazin & White (1987), 
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This fit is accurate to within 4% accuracy, for a plasma 
with solar metallicity, within 10 5 < T < 10 8 K. For 
10 s <T< 10 9 K, it underestimates cooling by a factor of 
order unity (compared to the exact cooling function , as 
in, e.g., Schmutzler & Tscharnuter, 1993) , and therefore 
is a conservative fit to use, as far as the effect of cooling is 
concerned. 

For non-steady flows, we adopt the formalism of White 
& Sarazin (1987), where the mass deposition rate, p, is 
characterised by a 'gas-loss efficiency' parameter q. One 
writes p = q(p/t coo i) where t coo i is the local isobaric cool- 
ing rate (t coo i = 5fc^T 'pm p / pA). It has been found that 
q ~ 3 models can reproduce the observed variation of mass 
flux (m cx r) (Sarazin & Graney 1991). Fabian (1994) has 
noted that these models of White & Sarazin (1987) yield 



good approximations to the emission weighted mean tem- 
perature and density profiles for cooling flow clusters. We 
also note that Rizza et al (2000) have used the steady 
flow models of White & Sarazin (1987) to simulate cooling 
flows. 

We first discuss cooling flows with m ^constant. With 
q = 0, one can eliminate the density from Eq. [ll] to get 
two differential equations: 
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These equations have singularities at the sonic radius r s 
where 59 s = 3u s . A necessary condition of singularity is 
that the numerators of Eq. vanish at the sonic radius. 
Therefore (Mathews & Bregman 1978) 
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We have used two different dark matter profiles for the 
cluster. The first model ( Model A) has been discussed 
earlier in the literature in the context of cooling flows in 
cluster (White & Sarazin 1987; Wise & Sarazin 1993) with 
a density profile, 
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Here p a = 1.8 x 10 gm cm and r core = 250 kpc de- 
scribe the profile of the cluster mass, and p o g = 4.1 x 10 -22 
gm cm -3 and r c g — 1.69 kpc describe the profile of the 
galaxy in the centre of the cluster. 

Model B does not have the galaxy in the center, and so 
it is described simply by p = po/[l + (r /r core ) 2 }. 

With these assumptions, the solutions for steady cool- 
ing flows, to = constant) are fully characterized by (1) 
the inflow rate, 777, and (2) the temperature of the gas T s 
at the sonic radius r s . Obviously, the cooling flow solu- 
tions are only valid within the cooling radius r coo i where 
tcooi = t age - We assume a value of t age = 10 Gyr for 
all models. We assume that outside the cooling radius, 
gas obeys quasi- hydrostatic equilibrium (Sarazin 1986). 
Although this means matching the cooling flow solutions 
with nonzero u to u = solutions outside, in reality the 
velocity of gas at the cooling radius is very small (for a 
777 = 100 M Q yr _1 with p - 10~ 26 gm cm" 3 , at r = 250 
kpc implies a velocity of 30 km s^ 1 ), which is close to 
the limit of turbulence in the cluster gas (Jaffe 1980), and 
smaller than the sound velocity (~ 1.5 x 10 3 (T/10 8 K) 1 / 2 
km s _1 ). The velocity of the flow at the cooling radius 
is, therefore, for all practical purposes, sufficiently small 
to be matched to the solution of hydrostatic equilibrium 
outside. (In this approach, we avoid the time consuming 
search for the critical value of 777 for which the flow so- 
lutions behave isothermally at r — > 00 (sec Sulkanen et 
al. 1993).) 

As in the usual assumptions for the interpretation of SZ 
effect, we assume that the gas outside the cooling radius 
is isothermal, with a constant temperature profile. The 
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density, therefore, obeys p oc [1 + (r/r core ) 2 ]~ 313 / 2 , where 
= fj,m p a 2 /kBT iso , and T iso is the temperature of the gas 
at and outside the cooling radius. 

For models with non-zero q (Model C: has the same mass 
profiles as Model A) , the solutions are characterized by T s 
and the value of m at the cooling radius, m coo i. Since a 
fraction of mass drops out of the flow in this case, the in- 
flow velocity need not rise fast and so it is possible to find 
completely subsonic solutions. 

3.2. Cooling flow solutions 

We numerically solved the flow equations for the pa- 
rameters listed in Table 1. The density, temperature and 
pressure profiles for three cases are presented in Figures 
1, 2 & 3 . We mark the position of r coo i in each case, 
and we mark r s for the cases of transonic flows (when 
m = constant). Beyond r coo i we match a hydrostatic solu- 
tion, as explained above, for the respective potentials. We 
also present, for comparison, the behaviour if the solutions 
outside r coo i arc assumed to extend inwards (that is, if no 
cooling flow is assumed). We will postpone the discussion 
on the effect of these profiles on the SZ decrement to a 
later section, and only discuss the qualitative aspects of 
the solutions here. 

The solution Al is similar to that presented by Wise 
& Sarazin (1993) (their Figure 1; although they chose to 
characterise the solutions by the temperature at r coo i and 
not T s as we have done here). It is also similar (qualita- 
tively) to the solution for A262 presented by Sulkanen et 
al. (1989). As the latter authors have noted, the effect of 
having a galactic potential in the center is to have a flat- 
ter temperature profile for r > r s , than in the case of no 
galactic potential. This aspect is clearly seen while com- 
paring our solutions with and without galactic potentials 
in the center. Our calculations for the case without the 
central galaxy are admittedly flawed in the very inner re- 
gions where the gas density is larger than the dark matter 
density, which results in an incorrect determination of the 
graviational potential in the inner region. However, this 
happens only inside a region ~ 25 kpc from the centre, and 
should not influence our final results by a large extent. 

A word of explanation for the pressure profiles is in order 
here. Naively speaking, it would appear that the pressure 
profile inside the cooling radius should have lower values 
than the corresponding case of hydrostatic equilibrium. 
The fact that it is not always so has been noted in the lit- 
erature (e.g., Soker & Sarazin 1988, Fig 1 of Sulkanen et 
al. 1989 ). The reason for the pressure bump just outside 
of the sonic radius is that the flow in this inner region is 
not pressure driven, but rather by gravity (see also Soker 
& Sarazin 1988). This is why the bump in the profile de- 
pends on the presence and absence of the galaxy in the 
centre. And this profile leads to the curious result that 
the presence of cooling flow can lead to the overestimation 
of the Hubble constant as discussed in the next section 

The model with mass deposition (CI) is shown in Figure 
3. The local mass flux is found to be almost proportional 
to the radius, consistent with various observations (Fabian, 
Nulsen & Canizares 1984; Thomas et al. 1987), and, there- 
fore, is probably a realistic model for cooling flow clusters. 
In this case the temperature drops gradually all the way 
through, since the velocity does not rise too fast. The de- 



posited mass is assumed to impart negligible pressure and 
the pressure refers only to that of gas taking part in the 
flow. 



4. DETERMINATION OF HUBBLE CONSTANT 

In this section we discuss the SZ and X-ray profiles of 
clusters with cooling flows. We compare these with profiles 
from cluster having gas in hydrostatic equilibrium, and 
comment on the reliability of measuring Hubble constant. 
The effect of cooling flow and the subsequent increased 
Brcmsstrahlung emmison is seen in the sudden increase in 
the X-ray flux in the innermost region of the cluster (Fig- 
ure 5). The signature of the cooling flow is seen in the form 
of the central spike in the X-ray profile. The X-ray profile 
is only affected slightly by the drop in temperature and 
it is the dependence on the gas density that holds . The 
temperature dependance becomes important only near the 
sonic point. Outside r coo ;, the X-ray profile is the same as 
that in the hydrostatic cases. 

The SZ distortion is proportional to the line of sight 
integral of the pressure, and the sudden increase of the 
gas density inside the cooling radius is moderated by the 
decrease of the gas temperature. As a result there is a 
gradual increase in the gas pressure. Near the sonic point 
the temperature drops drastically by orders of magnitude, 
and results in sudden decrease in pressure. However, since 
this change in pressure becomes acute only within w 5% 
of the core radius, it contributes negligibly to the line-of- 
sight integral of the gas pressure, and leads to an increase 
in the SZ distortion inside the cooling radius for all mod- 
els considered (see Figure 4). Like the X-ray profiles, the 
SZ profiles outside r coo i is the same as that for the corre- 
sponding hydrosstatic cases. 

The SZ profiles have been calculated in the Rayleigh- 
Jeans limit (x <C 1) where j v {x) of Eqn. 4 goes to —2. In 
general, however, the profiles should be calculated using 
Eqn 2. Our results below are independent of the observa- 
tional frequency, since the profiles at different frequencies 
have similar shapes, with the amplitude of the SZ distor- 
tion scaled either up or down. 

Once both profiles are known, one can determine the de- 
viation in the value of the Hubble constant using Eqns(14) 
& (16). The deviation from the idealistic case can be 
parametrised as 

r Tcore,true Ho.est /oo\ 

JH = ■ = 77 (22) 

'core, est, /«i *l0,true 

The above formula has been used to determine the devi- 
ation of the estimated value of H from the actual value, 
for models listed in Table- 1. The effect of cooling flow 
on the determination of the cosmological parameters are 
summarised in Table-2. 

To begin with, one has to get best fitted values for r core 
(or C ) and from different profiles. Since, the estimation 
of the Hubble constant depends on the determination of 
these parameters from the profiles, we look at this issue 
in more detail. We must keep in mind that the best fit- 
ted value of r core (or 6 C ) and depends on whether one 
decides to fit the X-ray or the SZ profiles, and the choice 
can lead to significant differences in the estimated value of 
H . One of the reasons for the strong dependance on the 
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nature of the profile can be the non-isothcrmality of the 
cluster gas. Recent observations indicate that intraclus- 
ter gas has a temperature structure, see Markevitch et al 
1998. This is because the y-parameter depends on the in- 
tegral over T e , while emissivity of thermal Bremsstrahlung 
depends on \JT e . The dependance of the Gaunt factor on 
T e isindirect and weak. Yoshikawa et al 1998 have shown 
that gas temperature drop in the central regions (their Fig. 
3), should increase both r core and (3 fitted to y(9), and to 
alesser extent to Sx{6), as compared to those compared 
to n e (r). This discrepancy increases at higher rcdshifts. 
However, in their case, there is little change in the gas 
density profile, dumpiness can also give rise to different 
fits, resulting in an overestimation of the Hubble constant 
(Inagaki etal 1995). 

There are two other important points that have to be 
kept in mind while fitting the profiles. First, we must re- 
member that we are trying to fit a cluster having a finite 
profile with the formulae (Eqns 12 & 13) for isothermal 
(3 profiles which is derived assuming the cluster to be of 
infinite extent. This can, by itself, lead to an overestima- 
tion of H (Inagaki et al 1995). Thus to have a good fit 
one must choose a segment of the profile such that, within 
that segment, the profiles (SZ or X-ray) for a finite cluster 
do not differ much from those of a hypothetical cluster of 
infinite extent. We found that SZ and X-ray profiles of 
clusters start differing from those of infinite size at radii 
greater than 1.5 times the core radius. Hence, we have 
restricted our fitting to radii within 1.5r core . 

Next, one must also be careful to exclude the region 
close to the sonic point, so that the X-ray spike is ex- 
cluded from the fit. Also, the central portion in the SZ 
profile should be avoided as its inclusion can give an ap- 
parent central decrement less than its neighbouring points 
(see Schlickeiser 1991). We have fitted the SZ and X-ray 
profiles varying the inner cutoff radius and the results for 
a representative solution for each class of model are tab- 
ulated in Tables 3 & 4. Thus, all fittings were done for 
profiles extending from r = r m i n to r = 1.5r core . 

As can be seen from Table-2, cooling flows can lead to an 
overestimation of the Hubble constant. However, we must 
emphasise, that it may not be possible to a priori estimate 
the amount of bias introduced in the measurement of the 
Hubble constant due to cooling flows. There is no simple 
correlation between the amount of cooling (i.e m) and the 
change in the estimated H from the actual value. The 
total change depends not only on tin, but also on the po- 
sition of cooling radius, sonic radius, temperature at the 
sonic point and the isothermal temperature characterising 
the hydrostatic cases, with which comparisons are made. 
Specifically, the fitted values of r core (8 c ) and (3 for cooling 
flow models differ from hydrostatic models according to 
shape of underlying profiles, which is marked by two im- 
portant features, firstly the central excess of X-ray flux (or 
excess decrement of SZ flux) , and secondly the the devia- 
tion from the smooth hydrostatic profile inside r coo i, the 
amount of overestimation mainly depends on these factors. 
For models with a central galaxy potential, there is always 
an over-estimation of H Q , which is greater than the models 
without the central galaxy. 

For the realistic cases of models CI and C2, where we 
have a variable m with r inside the cooling radius, the devi- 



ation of estimated Hubble constant from its actual value is 
almost the same. They are also greater than that of mod- 
els A and B, having similar mass flow rates. This may be 
due to the fact that the maximum deviation in pressure 
from the hydrostatic cases is more in non steady cases, 
than in steady flows. Also non-steady cases are marked 
by the absence of the sonic radii and the subsequent drop 
in temperature. 

We note that although the different choice of fitting may 
change the absolute determination of cosmological param- 
eters, the trend i.e deviation from the correct values re- 
mains more or less unaffected. It is interesting to note 
that for B-type model (CI), which include mass deposi- 
tion in cooling flows, the deviations decrease as one ex- 
cludes a greater part of the cooling flow region (Table 4& 
5). The other models show an increase instead. Here, we 
remind ourselves that models with mass deposition i.e C- 
type models are more realistic (Fabian 1994). It is possible 
that the unsually high value of deviation (Tabic 4) and 
the counter-intuitive trend of increasing deviation with 
decreasing portion of cooling flow region used for fitting 
(Table 4 & 5), arise because of the unrealistic modelling of 
cooling flows. If we take the model Bl as a realistic one, 
then Table 4 & 5 show that to obtain a value of the Hubble 
constant within an accuracy of <~ 10%, one should have 
r m in ~ 0.8r coo; . In most cases, r coo i < r core (Fabian et 
al. 1984). However, as r coo ; cannot be determined without 
actually detecting a cooling flow in a cluster, we suggest 
that a significant portion of the profile within r core should 
be excluded as a precaution. The SZ and the X-ray profiles 
for the different models are shown in Figures 4 & 5. 

5. DISCUSSIONS AND CONCLUSION 

Our work on the effect of temperature structure of clus- 
ters and its effect on SZ decrement differs from other previ- 
ous work of this nature in following way: this work takes 
into account the change in density profile as well as the 
temperature profile since both becomes important in the 
central region of the cluster. Also, previous authors have 
looked at the issue of non-isothcrmality of a cluster at 
radii greater than the core radius of the cluster, whereas 
we look at temperature change at regions inside the core 
radius. For them the density profile can still be well ap- 
proximated by a (3 profile, whereas for cooling flow solu- 
tions, density profile is vastly different. Further, they have 
neglected radiative cooling in their work. We for the first 
time look at SZ effect in presence of radiative cooling, by 
first solving the cooling flow equations for reasonable and 
physical solutions. 

In summary, we find that the presence of a cooling flow 
in a cluster can lead to an overestimation of the Hubble 
constant determined from the Sunyaev-Zel'dovich decre- 
ment. We have used different models of cooling flows, with 
and without mass deposition, and found the deviation in 
the estimated value of the Hubble constant in the case of 
a cooling flow from that of hydrostatic equilibrium. We 
have used the usual procedure of fitting the SZ and X-ray 
profile with a [3 profile to get an estimated value of r core , 
and then compared with that for the case of gas in hydro- 
static equilibrium in order to estimate the deviation in the 
Hubble constant. For the more realistic models with mass 
deposition (varying rh with radius), we found that the de- 
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viation decreases with the exclusion of greater portions of 
the cooling flow region. Quantitatively, we found that for 
the deviation to be less than ~ 10%, one should exclude a 
portion of the profile upto <~ 0.8r coo /. Since r coo ; is diffi- 
cult to estimate without actually detecting a cooling flow, 
we have suggested that a significant portion of the profile 
inside r core should be excluded, to be safe. 

There can be another important implication of the ef- 
fect of cooling flows. With the upcoming satellite mis- 
sions (MAP & Planck), we have come to the point where 
there are efforts to constrain f2 with surveys of blank SZ 
fields (Bartlett et al 1998; Bartlett 2000), ultimately giv- 
ing rise to SZ-selected catalog of clusters (Aghanim et al 
1997). This method relics on estimating the number of 
SZ sources brighter than a given threshold flux (Barbosa 
et al 1996). The point to be noted is that since these 
surveys are essentially flux limited in nature, the valid- 
ity of the analysis to determine f2o depends crucially on 
the one-to-one association of flux-limits to corresponding 
mass limits of clusters. From our analysis above, it seems 
that it may not be possible to associate a unique cluster 
mass to a given SZ distortion, given the uncertainty due 
to the presence of cooling flows. This might lead to con- 
tamination in SZ cluster catalogs and the inference of flo. 
Recently, attempts have been made to constrain Q from 
variance measurement of brightness temperature in blank 
fields (Subrahmanyan et al 1998) and comparing them to 
simulated fields (Majumdar & Subrahmanyan, 2000), of 
cumulative SZ distortions from a cosmolgical distribution 
of clusters. These results may also be systematically af- 
fected due to the presence of clusters having cooling flows. 

The estimations made in this paper strictly applies to 
cases where the image of the SZ effect is directly obtained 
by single dish observations. For interferometric observa- 
tions, the interferometer samples the fourier transform of 
the sky brightness rather than the direct image of the sky. 
The fourier conjugate variables to the right ascension and 
declination form the u — v plane in the fourier domain. 
Due to spatial filtering by an interferometer, it is neces- 
sary that models be fitted directly to the data in the u — v 
plane, rather that to the image after deconvolution. We do 
not forsee drastic change from our inferences in such cases 
since the result mainly depends on the deviation of the SZ 
and X-ray profile in case of a cooling flow from those in hy- 
drostatic equilibrium. This, however, should be looked in 
greater detail in future. We also note that with the grow- 
ing number of high quality images of the SZ effect with 
interferometers, which have greater resolution than single 
dish antennas, the shape parameters of the clusters can be 
directly determined from the SZ dataset rather than from 
an X-ray image (Grego et al. 2000). 

Finally, we would like to add, though the calculations 
presented in this paper were done using the dark matter 
profile (Eqn. 21), which is "commonly used" for calculat- 
ing cooling flow solutions, it is inconsistent with the dark 
matter profile (Navarro et al 1997) found in numerical sim- 
ulations. (For a comparison of mass and gas distribution 
in clusters having cooling flows with different dark matter 
profiles, see Waxman & Miralda-Escude', 1995). More- 
over, we have neglected the self gravity of the gas. Suto 
et al (1998) have calculated the effect of including the self 
gravity of the gas in determining the gas density profile. To 



make strong conclusions about the effect of cooling flows 
in the determination of the Hubble constant, one should 
take both the points mentioned above into account. 

We are grateful to Joseph Silk for his comments on 
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like to thank the Raman Research Institute, Bangalore, 
for hospitality and for providing computational facilities. 
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Table 1 - Parameters for cooling flow solutions 



Solution 


Mass Model 


q 


m{r cool ) 


T s 




f cool 


T 








(M yr-*) 


(K) 


(kpc) 


(kpc) 


(K) 


Al 


Model A 





100 


6.5 x 10 b 


0.688 


127.5 


1.2 x 10 8 


A2 


Model A 





200 


6.5 x 10 b 


0.462 


96.1 


7.7 x 10 B 


A3 


Model A 





300 


6.5 x 10 b 


0.712 


132.2 


7.7 x 10 8 


Bl 


Model B 





100 


4.0 x 10 b 


0.688 


85.7 


1.14 x 10'' 


B2 


Model B 





200 


6.5 x 10 b 


0.462 


89.6 


1.9 x 10'' 


B3 


Model B 





300 


6.5 x 10 b 


0.712 


110.3 


1.9 x 10 a 


CI 


Model A 


3 


200 






111.6 


1.1 x 10 s 


C2 


Model A 


3 


300 






132.2 


1.1 x 10 s 



Table 2 - Effect on central decrement and H for r m in = O.lr, 



Solution Type 










(% change) 


Al 


1.91 


35 


A2 


1.18 


11.5 


A3 


2.6 


25 


Bl 


1.36 


12.0 


B2 


1.19 


9.0 


B3 


1.13 


8.5 


CI 


2.6 


14.0 


C2 


2.7 


11.3 



Table 3 - Fitting of SZ profile and deviation of Ho for fa = 200m,Q/yr 



Solution Type 


Tmia — 0.2r coo / 


T'raia — 0.5r coo / 


T'min — 0.8r coo / 




H/ Htrue 


H/ Htrue 


H/ Htrue 


A2 


1.57 


1.74 


2.18 


B2 


1.44 


1.58 


2.06 


CI 


2.20 


1.18 


1.07 



Table 4 - Fitting of X-ray profile and deviation of H a for fa = 2OOm0 / yr 



Solution Type 


Train — 0.5r coo / 


Train — 0.8r coo / 


Train — 0.9r coo / 


Train — 0.95r coo / 




H/ Htrue 


H 1 Htrue 


H 1 Htrue 


Hf Htrue 


A2 




4.7 


2.7 


1.7 


B2 




4.9 


2.24 


1.6 


CI 


1.69 


1.12 


1.02 


w 1.0 
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-22 




Fig. 1. — Cooling flow solution A2. The upper left panel shows the dark matter density profile (dotted line), the gas 
density profile for the cooling flow (solid line) and the corresponding case of gas in hydrostatic equilibrium. The lower left 
panel shows the temperature profiles for the same cases. The position of r s and r coo i are shown. The upper right panel 
shows the pressure profiles and the lower right panel plots the Mach number of the cooling flow gas as a function of the 
radius. 




Fig. 2.- 



Same as figure 1, for cooling flow solution B2. 
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Fig. 3. — Same as figure 1, for cooling flow solution CI. 



solution type A2 



solution type B2 




0.081 
0.0805 
0.08 
; 0.0795 
' 0.079 
0.0785 
0.078 
0.0775 



solution type C2 




Fig. 4. — Unconvolved SZ profiles, for the solution types A2 (upper left), B2 (upper right), CI (lower left) and C2 
(lower right), for cooling flow (solid line) and thecorresponding cases of gas in hydrostatic equilibrium (dashed line). 
Tcmb ^ s Pitted against . Individual plots have been magnified to highlight the differences between cooling flow and 
hydrostatic cases. 
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Fig. 5. — Unconvolved X-ray profiles for the solution types A2 (upper left), B2 (upper right), CI (lower left) and C2 
(lower right), for cooling flow (solid line) and the corresponding cases of gas in hydrostatic equilibrium (dashed line). 
X-ray surface brightness in units of erg s~ x cmr 2 ster~ x is plotted against — r —. Individual plots have been magnified to 

T core 

highlight the differences between cooling flow and hydrostatic cases. 



